Fractional vortices on grain boundaries — the case for broken time reversal symmetry 

in high temperature superconductors 
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We discuss the problem of broken time reversal symmetry near grain boundaries in a d-wave 
superconductor based on a Ginzburg-Landau theory. It is shown that such a state can lead to 
fractional vortices on the grain boundary. Both analytical and numerical results show the structure 
of this type of state. 



I. INTRODUCTION 

During the last few years the understanding of the mi- 
croscopic properties of the high-temperature supercon- 
ductors (HTSC) has gradually improved While for a 
long time studies have focused on the exotic normal state 
properties, recently the interest turned more towards the 
superconducting phase, in particular, the symmetry of 
the order parameter. For a tetragonal system the list of 
possible order parameter symmetries is rather long 
However, the recent debate has essentially concentrated 
only on two symmetries of the Cooper pair wavefunction 
||. One is due to "s-wave" pairing, the most symmetric 
pairing channel. The other is "c?-wave" pairing, where 
the pair wavefunction (ip(k) oc cosk x — cosk y ) changes 
sign under 90°-rotations in the basal plane of the tetrag- 
onal crystal lattice. As a consequence, the latter wave 
function has nodes along the [110]-direction. A possible 
alternative to the standard s-wave was presented with the 
"extended s-wave" pairing state (V>(k) oc cos k x + cos k v ) 
which has also nodes in the first Brillouin zone, but is 
completely symmetric under all operations of the tetrag- 
onal point group Q. For the orthorhombically dis- 
torted system, the s- and d-wave channels are not dis- 
tinguished by symmetry. Nevertheless, we expect that 
basic properties of the pair wave function such as the 
existence of sign changes and nodes are retained if they 
were present in the tetragonal case Q . 

A variety of experiments have been performed in or- 
der to distinguish among the order parameter symme- 
tries. One class of experiments considers the proper- 
ties of the quasi particle excitations in the supercon- 
ducting state. The existence of nodes in the pair wave- 
function implies that there are also nodes in the exci- 
tation gap. Low-lying excitations at the nodes modify 
the low-temperature behavior of certain thermodynamic 



properties compared with that of a superconductor which 
opens a complete gap. The clearest sign of such an effect 
was observed for the London penetration depth which 
behaves as A(T) — A(0) oc T in contrast to the conven- 
tional exponential law, oc exp(— A/fceT) ||. This re- 
sult strongly suggests that there are nodes in the gap 
and the pair wavefunctions, which are compatible with 
both extended s-wave and d-wave as well as with a very 
anisotropic s-wave state. 

Another class of experiments aims at the direct obser- 
vation of the intrinsic phase structure, the sign changes 
of the pair wavefunction. The Josephson effect as a phase 
coherent coupling of the order parameters of two super- 
conductors provides the natural means for this purpose 
Hf§. Arrangements connecting YBCO single crystals at 
two perpendicular surfaces to a standard s-wave super- 
conductor to form a loop for a SQUID have been used to 
detect a phase difference between the x- and y-direction 
of the pair wavefunction . The experiments observe 
with good precision a phase difference of tt compatible 
with the d x 2_ y 2--w&ve order parameter. 

The intrinsic tt phase shift in this configuration leads to 
frustration effects which manifest themselves in the form 
of a spontaneous supercurrent flowing around the loop. 
The supercurrent generates a flux $ = ±<E>o /2 where 
$0 = hc/2e is the standard flux quantum. This property 
has recently been detected and the flux was measured 
with very high accuracy iflTf . 

On the other hand, several other experiments based 
on the Josephson effect seem at present to contradict 
the presence of a d-wave order parameter. Chaudhari 
and Lin analyzed the Josephson current through a grain 
boundary in the basal plane with a special geometry giv- 
ing a basal plane contact between two segments of a 
YBCO film |Q. They demonstrated that various prop- 
erties might support an order parameter with s-wave 
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rather than simple d-wave symmetry. The interpretation 
of this experiment, however, has recently been contested 
by Millis |l|]. In contrast, Sun and coworkers investi- 
gated Josephson tunneling between a standard s-wave 
superconductor (Pb) and YBCO, where the tunneling 
direction is the c-axis of YBCO [ [L4| |. These data so far 
could not be explained consistently within the picture 
of pure d-wave superconductivity. Therefore, the sim- 
ple d-wave scenario may not be sufficient for a complete 
understanding of all experiments introduced here |15| . 

Indeed a recent experiment by Kirtley and cowork- 
ers suggests that the situation is more complicated than 
might be naively expected for a d-wave superconductor 
pa . Their experimental arrangement consists of two seg- 
ments of c-axis textured YBCO films where one is a tri- 
angular inclusion within in the other. The basal plane 
crystalline axes are misaligned with one another. The 
boundary of the triangle acts as junction between the 
two segments. We will show in Section II that if YBCO 
were a d-wave superconductor we would expect vortices 
to appear spontaneously at two of the three corners of 
the triangle each containing a flux of ±<I>o/2. The ex- 
periment does find spontaneous vortices at corners, but 
these vortices have fluxes different from n$o/2 (n: inte- 
ger). In addition, flux appears at all three corners and 
occasionally also on an edge of the triangle. We will 
argue in Section II that this can be explained by a super- 
conducting state which violates time-reversal symmetry 
T. Therefore, the simple picture of a single component 
d x 2 _ y 2 -wave order parameter might not apply here. 

T-violation is not uncommon in the field of uncon- 
ventional superconductivity. A large number of super- 
conducting states classified by symmetry indeed break 
time-reversal symmetry . In the complex supercon- 
ducting phase diagrams of the heavy fermion compounds, 
UPt 3 and U!- x Th x Be 13 (0.02 < x < 0.045) states ap- 
pear which probably break time-reversal symmetry. It 
was shown theoretically that such superconducting states 
can generate spontaneous supercurrents and magnetic 
field distributions in the vicinity of lattice defects and 
surfaces {H^. In both compounds the occurrence of such 
local fieldsin connection with the superconducting phase 
transition has been detected by means of muon spin ro- 
tation (/LtSR) measurements |L9|. For both compounds, 
consistent phenomcnological theories for this effect have 
been formulated Jl7| . 

In the field of HTSC, various theories and mecha- 
nisms leading to T- violating superconducting states have 
been proposed. The effective two-dimensionality of the 
cuprates may serve as a basis for par ticles with fractional 
statistics, the so-called anyons EQ]. Laughlin showed 
that the resulting superconducting state has a composite 
order parameter of the form d x 2_ y 2 + ied xy which ob- 
viously breaks time-reversal symmetry Alternative 
mechanisms can lead to a T-violating states with the 
symmetry s + id x 2_ y 2 22 1. At present there is no in- 
dication beyond any doubt that such states are realized 
in the HTSC |p3fl . On the contrary, recent experiments 
demonstrate that at least at the onset of superconduc- 
tivity (T « T c ) the critical behavior of the London pene- 



tration depth is that of a single component order param- 
eter belonging to the universality class of the XY-spin 
model (24]. Only below an additional superconducting 
transition at lower temperature could the composite T- 
violating order parameter appear. No signs of such an 
additional phase transition have yet been observed in 
the thermodynamic properties. In addition, it should 
be noted that each of the T- violating states mentioned 
above lacks gap nodes. This fact would also lead to in- 
consistency with low-temperature measurements of the 
London penetration depth which demonstrate the pres- 
ence of nodes, as mentioned above ||. 

In this paper we will show that there is no conflict 
between the interpretation of the experiment by Kirtley 
and coworkers jl6| which could indicate T-violation and 
the other experiments which obviously rule out the exis- 
tence of such a state |B5[ . We argue that the latter ex- 
periments address bulkproperties, while the former one 
considers effects in connection with interfaces and grain 
boundaries. The seeming conflict is resolved when we 
assume that T-violation occurs only locally in the im- 
mediate vicinity of an interface. The bulk, on the other 
hand, may only have a single component order param- 
eter, presumably with d-wave symmetry, but we cannot 
rule out other symmetries. As we will discuss below, the 
extension of the T-violating state towards the bulk is 
rather short, of the order of coherence length £. 



II. A FIRST INTERPRETATION OF THE 
EXPERIMENT 

Let us now examine the properties of an arra ngem ent 
similar to the one used by Kirtley and coworkers |16| . As 
illustrated in Fig. la, it is a superconducting film of tri- 
angular shape as an inclusion in another superconducting 
film, both of the same material. The crystal symmetry is 
tetragonal (for simplicity we neglect here the orthorhom- 
bic distortion present in many HTSC) and the film is c- 
axis textured. The basal plane axes (of the inclusion and 
the surrounding) are misoriented with each other. The 
interfaces (the edges of the triangle, each of length L) are 
weak links between the inner and the outer film. For sim- 
plicity we will treat them as Josephson contacts so that 
the standard sinusoidal current-phase relation applies. 



A. Pure d-wave symmetry 

Let us analyze the properties of this arrangement un- 
der the assumption that the superconductor here is a 
d-wave superconductor with an order parameter symme- 
try as the pair wavefunction ^d(k) = cosfc x — cosk y . 
This means that we should carefully consider the intrin- 
sic phase structure of the order parameter when deriving 
the Josephson current-phase relation. The phase differ- 
ence between the positive and negative lobes of the pair 
wavefunction is 7r. If dominant lobes of the same sign face 
each other at an interface, the corresponding Josephson 
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current-phase relation has the standard form and the in- 
terface energy is minimized by a vanishing difference the 
order parameter phases (0-junction). However, if the fac- 
ing lobes have opposite sign, an additional phase 7r enters 
and the energy is minimized by a phase difference of it 
(7r-junction) H. We have 

Ej((p) = -^^-cos{(p-a) (1) 

where a = for a 0-junction and n for a 7r-junction, and 
ip is the phase difference through the interface. In Fig. 
la we assume that the edge segments 1-2 and 3-1 can be 
labeled as a 7r-junction and segment 2-3 as a 0-junction. 
This definition is not unique. A redefinition of the order 
parameter phase in one of the two superconductors (jfi — > 
<f> + 7r) would reverse this labeling. 
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FIG. 1. Triangular grain boundary in a d-wave supercon- 
ductor: (a) The edges act as Josephson junctions with phase 
shifts, or 7r; (b) The phase p tries to be pinned at phase shift 
values (indicated by dashed lines) and change in an antikink 
(kink) of width Aj at the corners 2 and 3. There are several 
possible solutions for ip due to the 2-7r-periodicity of the phase. 
The solution with an anti-kink and a kink leads to an overall 
phase winding around the triangle. An equivalent solution 
with two kinks (also shown) would give 2ir winding. 

We now map all segments of the interface onto a one- 
dimensional axis with periodic boundary conditions for 
the coordinate x as shown in Fig. lb (x + 3L = x). Here 
we can study the spatial variation of ip along x by using 
the Sine-Gordon equation 

d%<p = A7 2 sin((^-a(i)), (2) 

where both the Josephson penetration depth Xj = 
{<f>oc/ 'S^dlc) 1 ' 12 (d; magnetic width of the interface) and 



the intrinsic phase shift a are assumed to be constant 
within each segment (see Fig. lb) pf|. For good junc- 
tions (Aj "C L), (p tends to be pinned to the value a 
in each segment of the interface, but has to change at 
the boundaries where a is discontinuous. The solution 
of Eq.(2) shows kinks at these boundaries with an exten- 
sion of A j (Fig. lb). Note that the two kinks at 2 and 3 
can be either "kink" and "anti-kink" or both "kinks" as 
a consequence of the periodicity of Eq. (2). Other types 
of kink solutions arc energetically more expensive. 

The spatial variation of ip induces a local magnetic 
flux density on the interface given by the expression 
4>(x) — &od x ip(x)/2ir Therefore, each kink corresponds 
to a local flux line or vortex with a magnetic flux <£> = 
^oi'Pr — l Pi)/2ir, where tp r (i) denotes the values of tp far 
enough to the right (left) of the kink such that d x ip is 
essentially zero. Both kinks in Fig. lb indicate vortices 
with $ = ±$o/2. Due to the periodicity of x we find that 
ip(x + 3L) = ip(x) + 2irn (n: integer), which requires that 
the total flux integrated over the whole triangle interface 
be an integer multiple of $o- Of course, the triangle is 
surrounded by a superconductor whose single valued or- 
der parameter allows phase windings of 27m only. This 
"sum rule" implies that half-integer flux lines can only 
appear at two of the three corners. Because the flux on 
each corner can only vary by n$o ; at corner 2 and 3 there 
is always a flux line with a flux of at least $o/2. Larger 
fluxes could be stabilized by an external field. This result 
is equivalent to the one presented in Ref. 8 and 13. 

The comparison of our "experiment" with the one per- 
formed in reality shows that the simple picture we tried 
to draw here does not explain the measurement by Kirt- 
ley and coworkers Jl6| . They found fluxes at all corners, 
all of which are clearly smaller than $o/2- We call them 
fractional vortices. In all samples checked, the sum rule 
constraining the total flux on the boundary to an integer 
multiple of (f>o was satisfied with good accuracy. 

In the introduction we claimed that the existence of 
fractional vortices requires a superconducting phase with 
broken time-reversal symmetry. We give here a brief 
argument for this statement. Consider one of the cor- 
ners of the triangle (or a similar structure) with a vortex 
whose flux is $. Apply the time-reversal operation to 
this system. This reverses the flux ($ — > — $). If the su- 
perconductor is otherwise invariant under this operation 
(T-invariant), the difference between $ and — $ must be 
an integer multiple of $o as in every standard supercon- 
ductor: $ = — $ + n$ - Therefore, $ (= n$ /2) is 
an integer or half-integer quantum of $o as seen above. 
Consequently the observation of a vortex with a flux dif- 
ferent from those values can only mean that the super- 
conducting state is not invariant under the time-reversal 
operation pTy . 

B. Josephson effect for a T-violating interface 

A T-violating superconducting order parameter con- 
sists of at least two components (e.g., d x i_ y i + id xy or 
s + id x 2_ y 2). We therefore restrict ourselves to the case 
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of a two-component order parameter with a generic pair 
wave function ^(k) = 771^1 (k) + 772^2 (k). Here 7/1 and 772 
are the two complex order parameters with the symme- 
try properties of the corresponding pair wavefunctions. 
Time-reversal transforms the order parameter to its com- 
plex conjugate (r/j — > 77*). If time-reversal symmetry is 
conserved, then (771,772) is up to a common phase fac- 
tor equal to (771,772) or 771/772 = Vi/vi- Otherwise, the 
order parameter breaks time-reversal symmetry and the 
state is at least two-fold degenerate, since both (771,772) 
and (771, 772) have the same free energy. We will consider 
here the standard situation where 771/772 is imaginary: 

- <t>2 = ±tt/2 {rjj = |?7j|exp(i0j))- 

For an interface between two superconductors (A and 

B) both with order parameter (r)f^ B \ vt^) * ne Joseph- 
son current phase relation has the form 



J = J cij sin( 



(3) 



where J is the supercurrent density at a given point on 
the interface. There are four different combinations for 
the phase coherent coupling and J C ij denotes the coupling 
strength between the components 77^ on side B and with 
rjf on side A (rjj — \r^\ exp(i^)) with j = 1,2 and 
fi = A, B). The energy for a uniform interface of area 5* 
is given by 



E, = - 



£o5 
2ttc 



COS 



4) 



(4) 



For simplicity we restrict ourselves to the situation where 
<f>i — <j>2 = <f>f — <j>2 = tt/2 is fixed on both sides of the 
interface. (This phase difference can be different in the 
vicinity of the junction without changing the conclusion 
we will draw here for the simplified case.) The current 
density J and the interface energy density ej depend only 
on one phase difference through the interface, say ip — 



with 



J(<p) = J c sin(</3 — a) 



e j(f) = cos(<p - a) 



Jc — V (Jell + Jc22) 2 + (Jcl2 — Jc2\) 
Jcl2 — Jc21 



(5) 



(6) 



tan(a) 



Jell + J, 



1-22 



The phase shift a corresponds to tp minimizing the inter- 
face energy. In this sense the correct solution for a must 
be chosen in Eq. (6) ||. Obviously, a can assume any 
value and depends only on the relative magnitude of the 
different coupling components which parameterize here 
the interface properties. It is easy to follow the same 



consideration for the T-invariant combination of 771 and 
772 where we find that a is strictly either or n. 

Let us apply this result to the triangle studied above, 
assuming T-violating superconducting states. Each of 
the interface segments as a uniform junction is charac- 
terized by a phase shift a (0 < a < ir). Analyzing Eq. 
(2) for this situation we find that there are now kinks of 
ip at all three corners connecting the different values of a 
in each segment (Fig. 2). This leads to vortices at these 
corners whose fluxes are given by 



$ = $, 



at 



Q- 



2tt 



(7) 



with a r m as the values of a on the segments on right 
(left) of the corner. It is obvious that in this case the flux 
$ need not be $0 or ^o/2< Thus, these kinks correspond 
to what we described as fractional vortices above. Note, 
that also here the sum of all fluxes must add up to an 
integer multiple of $0 because of the periodic boundary 
condtions for (p. 



9 

71 







OC12 



1 




OC31 



«23 



X 



1 



FIG. 2. Phase ip for a triangle with T-violating interface 
states. The phase shifts are indicated by the dashed line. 
tp follows these phase shifts by creating (anti-)kinks at all 
three corners. The kink heights are different from a multiple 
of 7r in general and lead to fractional fluxes. 



III. THE INTERFACE STATE 

As mentioned earlier, the assumption of a T-violating 
bulk superconducting state is incompatible with a num- 
ber of experiments. It is, however, possible that T- 
violation occurs only in the vicinity of the interface (grain 
boundary). This is enough to generate fractional vortices 
at grain boundary corners. In this section we would like 
to discuss such an interface state based on a Ginzburg- 
Landau theory. 



A. Ginzburg-Landau theory 

We use a Ginzburg-Landau free energy functional of 
two complex order parameters 771 and 772 . The first order 
parameter shall be the dominant one corresponding to 
the d x 2_ y 2-w&~ve pairing state. For the second we choose 
the ctcy-wave pairing state. Thus the pairing state is a 
combination of the two components 
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V>(k) = Yl wkOO = Wi (k) + ^2 (k) (8) 

i=l,2 

with 

^i(k) = fc2_ fc 2 and ^a(k) = Mfc. (9) 

The free energy functional must be invariant under the 
symmetry transformations of the crystal group (in this 
case tetragonal, D^), time reversal and U(l)-gauge sym- 
metry. Keeping terms up to fourth order in the expansion 
with respect to order parameter, we have 

T = J d 2 x \Ti + T-a + Ei=i,2^i] (10) 

where 

Ti = a t \ m \ 2 + + Ki\U m f (11) 

^12 = + 6{vh? + cc.) + k 2 (V x A) 2 

+ £ ((iVh)*(JV/2) - (D yVl y(D xm ) + cc.) (12) 

i=l,2 

- ^ J^V?r]j A }Hinterface). (13) 

We work in units where D = V i27rA/$o (A: vec- 
tor potential). The coefficients are all phenomenological 
parameters which contain the relevant information of mi- 
croscopic origin. The second order coefficients depend on 
the temperature in the usual way (a 8 - = a' (T — T C i) where 
T c i is the bare transition temperature of rji). The inter- 
face free energy Ti includes the pair breaking effect in 
the first term and the Josephson junction energy in the 
second term (note: J12 = J2i). The superscripts B and 
A denote the different sides of the interface, or "outside" 
and "inside" if the interface is a closed curve. 

It is worth noting that the free energy expansion with- 
out interface term is cylindrical symmetric with z as the 
rotation axis, although the system has only tetragonal 
symmetry. Anisotropy enters via the interface term. The 
coefficients gi and Jjj depend in general on the angle be- 
tween crystal axes and the interface. Pair breaking for 
the d x 2_ y 2-w&ve (d xy -wELve) phase is most effective when 
the interface normal vector points along the [1,1] ([1,0])- 
direction, i.e. along the nodes of the pair wavefunction 
p7[ . On the other hand, pair breaking is weaker when 
the lobes of the pair wavefunctions point towards the 
interface. Consequently, the d x 2_ y 2-w&ve component of 
the order parameter is suppressed weakly where the d xy - 
wave component is affected most (and vice versa). Sim- 
ilarly, the interface coupling coefficients depend on the 
internal structure of the pair wavefunction. If we denote 
the angle between the interface normal vector and the 
crystalline x-axis by they are J,j = Jofi{$B)fj($A), 



where the functions have the generic angular structure 
/i($) = cos(2z?) and = sin(2i?) which indicate the 

internal phase structure of the two pair wavefunctions. 

The variation of T with respect to r\i and A yield the 
following Ginzburg-Landau equations: 

KiT) 2 T]i = a.^ + f3i\r)i\ 2 r)i + j\r) 3 -i\ 2 r)i 

+ 25 V 2 ^* - {-IfKB^-i (14) 

and 

k 2 V x (V x A) = n 2 j = J2 KiQmritVTH - |r/ 4 | 2 A) 

i=l,2 

+ KV x (zlmrnrfe) (15) 

The term Ti generates the boundary condition at the 
interface with normal vectors n A ' B , 

) i=i,2 

Kin B . (v _ ^ A)v B = g B v B -^-Y J ijV f (16) 
) J=l,2 

The imaginary part is related to the expression for the 
supercurrent j perpendicular to the interface and leads 
to the Josephson current (Eq. (3)) with 

Jaj = JMHl (17) 

where the order parameter values are taken at the inter- 
face. 

Before discussing the interface problem, we consider 
the properties of the order parameter in the bulk. We 
assume that T c \ > T C 2 such that for temperatures imme- 
diately below T c \ only the component 771 becomes finite, 
while r\2 remains zero. The instability condition for the 
occurrence of 772 is given by 

a 2 (T*) + (7- 21*1)177! (m 2 = 0, (18) 

defining the transition temperature T* < T c \. This sec- 
ond transition leads to a state where both 771 and 772 are 
finite. The relative phase, 6 = 4>\ — 02, depends on the 
sign of 8. For S < the combination is real, 8 = or ir 
and for 5 > 6 = ±7r/2. The latter state breaks time 
reversal symmetry. For the following we will assume that 
S > 0. However, the other parameters shall be chosen so 
that T* < in order to avoid the second transition as it 
is not observed in the experiment. We also require that 
7 > 26 so that the two order parameter components tend 
to suppress each other. 



B. Order parameter 

We consider now an infinitely extended interface. In 
this case the order parameter and the vector potential 
only depend on the coordinate perpendicular to the in- 
terface. With certain simplifications a qualitative discus- 
sion is possible as we showed in Ref. |25[ . We would like 
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to present here an analytical study and then substantiate 
the result by a complete numerical treatment. 

The basic concept leading to an unconventional inter- 
face state is the following. Pair breaking at the interface 
reduces the 771-component locally. It recovers, however, 
over a coherence length £ = yjK\/\ai\. It is easy to see 
from Eq. (18) that a local reduction of 771 leads to a local 
enhancement of T* (note: 7 > 28). Consequently, the 
772-component can appear at the interface at sufficiently 
low temperature (T < T'), but it decays exponentially 
towards the bulk. The extension of 772 diverges when T 
approaches the bulk T* . Therefore 772 does not possess 
the same length scale as 771 in general. Because 8 > 0, 
the combination of 771 and 772 is complex 9 7^ or tt. 

Let us now consider a simplified analytic solution of 
the Ginzburg-Landau equations for an interface whose 
normal vector is parallel to the ic-axis. We neglect the 
vector potential and the coupling between the two sides 
(Jij = 0). A qualitatively good view of the interface state 
is obtained in the limit where the length scales of the 
order parameters are very different, i.e. for £1 <C £2- At 
the interface, 771 behaves approximately like 



jyi(x) = ?7iotanh 



x 



6 



(19) 



as can be found by solving the Ginzburg-Landau equation 
with 772 = 0. The boundary condition at the interface 
determines xq and 7710 = y/\a*\/2/3i is the bulk value of 
the order parameter. Fixing rji (x) the equation for 772 
becomes 

K 2 d 2 xV2 = {a 2 +jT 1 2 1 (x))r l2 + 2Sr 1 2 1 (x)r 1 * 2 

+ P2\v2\ 2 m+g2mS{x) (20) 

Following our assumption about the coherence lengths, 
we approximate the spatial dependence of 771(2;) by a 8- 
function in this equation. It is easy to see that 772 is 
purely imaginary, 772 = iu(x) and satisfies the equation 

K 2 d 2 x u = a*u + fou 3 + (g 2 - (7 - 28) V l s^)8(x) (21) 

The factor s is of the order one. The solution is a hyper- 
bolic function. 



from the fact that we do not resolve here variations on 
length scales of £1. 

For finite Jij we can determine the phase shift a. In 
this symmetric formulation of the interface problem the 
phase shift a would be strictly or 7r even if the state 
is a complex combination of 771 and 772, i.e. it breaks 
time reversal symmetry. We find immediately that Eqs. 
(6) and (17) lead to a — or tt, because J\ 2 — J 2 \. 
A condition sufficient for a different from these trivial 
values would be different coefficients gt on both sides of 
the interface, which means for example that the crystal 
orientation on side A and B is different (as it is for a 
grain boundary). It can be shown that the violation of 
"parity" , the mirror symmetry due to reflection at the 
interface, is necessary (see M). 



C. Magnetic properties 

The interface state is an inhomogeneous T-violating 
superconducting state for temperatures below T'. It 
shows unusual magnetic features which originate from an 
orbital magnetic degree of freedom of the Cooper pairs. 
Remember that the two pairing components we consider 
both belong to the d-wave channel. If they are combined 
with a relative phase different from 6 = 0, 7r, then the 
pairing state has a component belonging to the spherical 
harmonic 5^,±2(k) oc (k x ± ik y ). This Cooper pair state 
has a finite orbital magnetic moment Mz parallel to the 
z-axis. The moment generates circular currents in the 
x-y-plane which cancel in the uniform superconducting 
phase. However, in the inhomogeneous region at the in- 
terface they can appear as finite currents, j = V x zM(r). 
These currents are included within our phenomenological 
description. 

Let us discuss the magnetic part of the Ginzburg- 
Landau equation, Eq. (14). Due to the translational in- 
variance parallel to the interface, no currents are allowed 
to flow through the interface in the energetically lowest 
state (otherwise, the Josephson energy would not be min- 
imized). Through the equation for the ^.-component, we 
find that this requires A x = 0. The ^-component then 
satisfies the following equation, 
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where £| = K 2 /a*. We use the boundary condition at 
the interface to find the shift xq (> 0), 
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Setting the right hand side equal one determines the criti- 
cal temperature T" for the occurrence of 772 . At this point 
xq = 00. Note that T" > T* and that 772 extends into the 
whole superconductor when we approach T — T* . In the 
range T* < T < T' it decays exponentially on the length 
£2- Our numerical results below will show that this ap- 
proximate solution describes the interface state well apart 
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where A~ 2 = K 2 t] 2 /k 2 is the London penetration depth. 
The right hand side denotes the current due to the mag- 
netic moment, which is proportional to ii(j]\ri 2 — rjirj 2 ). 
This current j y flows parallel to the interface. We do not 
discuss the influence of the vector potential on the order 
parameter here. Inserting the solution found above for 
the interface state, we find 
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This current distribution is odd under reflection through 
the interface. Starting at zero on the interface, it rises 
quickly and has a maximium at a distance of about £1. 
Then it changes sign and decays on the length scale £2- 
The magnetic field generated by this current has a narrow 
peak of width £1 on the interface followed by two wings 
of opposite sign. Within our approach we obtain 

B z = ^(u 2 (x) - sS(x)) (26) 

where we neglect the screening effects due to the sec- 
ond term on the left hand side of Eq.(24). Because the 
field in Eq.(26) leads to a finite magnetic flux, screening 
currents are induced which yield a compensating diamag- 
netic field with the length scale A. This screening effect 
can be rather small if the two contributions in Eq.(26) 
nearly cancel each other. The net magnetization vanishes 
exactly, because in the interior of a superconductor phase 
coherence allows for a net magnetic flux only if there is 
a winding of the order parameter phase. 
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FIG. 3. Numerical solution of the infinite interface state. 
Top: The modulus of the two order parameter compo- 
nents. The relative phase is constant 7r/2. Bottom: The 
z-component of the magnetic field (multiplied by k 2 ). In 
these units, £1 = 1 and A = 7 (k 2 = 49). The magnetic 
field exhibits variations on the lengthscales £i and £2 in the 
immediate vicinity of the interface, and is screened away from 
the interface on the lengthscale A. 



D. Numerical solution for the infinite interface state 

Finally we solve the complete set of coupled GL equa- 
tions for the uniform interface in order to show that our 
analytic treatment gives the correct qualitative behav- 
ior. We choose the coherence length £1 = y/K\/\a\\ as 
the unit length and the London penetration depth about 
10 times £1. For the coefficients we use values of order 
unity: a' = j3 2 = ft = K 2 = K x = 1, K = 0.1,7 = 1 
and 5 = 1/4. The transition temperature of the d xy - 
componcnt is T c2 = 0.4T c i. For the space coordinate we 
introduce a fine mesh and the interface is taken as single 
point. For simplicity we consider the symmetric situa- 
tion by representing the interface via a local suppression 
ofT c i: gi\vi\ 2 S(x). 

The solution for gi = 4, g 2 = is shown in Fig. 3. The 
shapes of the order parameters (top) are in good qualita- 
tive agreement with the analytic result. We see that the 
length scales £1 and £2 arc different and that the latter is 
larger than the former. The relative phase 9 between the 
two components is constant, 9 = tt/2. The magnetic field 
has a narrow peak with a width of the order £1 on the 
interface (bottom). Towards the bulk, the field changes 
sign and decays on the length A (= 7£i). The positive and 
negative parts of the field distribution cancel each other 
such that no net magnetization is present, as anticipated 
above. 



E. Numerical solution for the grain boundary corner 

In Sect. II we showed that at the border between 
two different grain boundary segments with T-violating 
states a vortex with fractional flux can appear. We would 
like to demonstrate this fact here by solving the complete 
GL problem for a system with a grain boundary which 
has a corner. For simplicity we use a right angle corner 
which matches well with our choice for the 2D-square 
lattice mesh. A variable grid has been introduced to en- 
hance the accuracy in the vicinity of the interface where 
the order parameter and magnetic field have the largest 
variations. A steepest descent (relaxation) method was 
used to minimize the GL free energy with open boundary 
conditions at the borders of the mesh. 

The coefficients of the GL free energy are the same 
in both grains (see caption of Fig. 5). However, the 
interface terms for the grain boundary are different for 
the two segments separated by the corner. The geometry 
of the grain boundary configuration is shown in Fig. 4. 
The A and B superconducting regions are separated by 
the grain boundary, represented as a thick black line. For 
the segment of the interface parallel to the xa~ (2M~)axis 
we choose J n = 0.1, J 12 = -(+)0.2, J 21 = -(+)0.2 and 
J22 = 0.1. In region A, gi = 2.5 and g 2 = 0. In region 
B, gi = 5 and g 2 — 0. As discussed above, the difference 
in these boundary conditions yields different phase shifts 
a in the two segments. For the parameters chosen, the 
phase shifts are equal in magnitude and opposite in sign 
in the two segments, as expected from Eqs. 6 and 7. 
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FIG. 4. Grain boundary configuration used in the calcula- 
tion of fractional flux. The grain boundary is indicated by 
the thick line. The A region and B region represent different 
crystal grains, characterized by their respective principal axis 
vectors as shown. The normal vectors at each interface are 
also shown. In the numerical calculation, the density of mesh 
points is greatest near the grain boundaries. 

In Fig. 5 we show the magnetic field distribution 
around the corner which has a pronounced peak indi- 
cating the position of the vortex. The flux is fractional 
of the size $ « —0.15 $o- Along the grain boundary the 
small field peak occurs which we found already for the 
uniform boundary case. This peak is distorted somewhat 
by the Josephson coupling of the order parameters across 
the interface. We also see a slight difference in the decay 
of the magnetic field towards the grains and along the 
grain boundaries. The grain boundaries are usually good 
contacts so that the magnetic flux at the corner is well 
defined. The sign and magnitude of the flux can be ma- 
nipulated by tuning gi and Jij. Setting all the Jij to 
eliminates the flux, of course. Switching the signs of J\i 
and J21 while leaving Jn and J22 unchanged reverses the 
sign of the magnetic flux, as expected from Eqs. 6 and 
7. Interestingly, by choosing the Jij and gi appropriately 
we can favor the formation of a domain wall along the 
grain boundary at one or both interfaces. The magnetic 
field distribution along the grain boundary in this case 
becomes antisymmetric with respect to reflection across 
the boundary. The parameters chosen for Fig. 5 do not 
lead to such a domain wall on either interface. 

The grain boundaries themselves are naturally good 
conduits for magnetic flux. In the case shown in Fig. 5, 
negative magnetic flux pours into the central dip in B z 
along the interface and drags down the side lobes (Fig. 3, 
bottom). In general, we expect the grain boundary to pin 
bulk vortices ($ = ±n$o) on top of the fractional vor- 
tices due to this tendency to grab magnetic flux. This is 
consistent with the experimental observations of Kirtley 
and coworkers [fL6|. 
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FIG. 5. The magnetic field distribution at this grain 
boundary corner. Top: contour plot. Lighter (darker) shades 
indicate more positive (negative) magnetic fields. Bottom: 
surface plot of — k B z . The parameters used are £1 = 1, 
k = 7, K x = Ki = ai = pi = fa = 7 = 1, <5 = 1/4, K = 0.1, 
and k = 7. For the segment of the interface parallel to the xa- 
(y a-) axis, Jn = J22 = 0.1, J12 = J12 = — (+)0.2. In region 
A, gi — 2.5 and Q2 = 0. In region B, gi — 5 and 172 = 0. The 
enclosed flux is $ « —0.15 $o- 



IV. CONCLUSIONS 

We have studied the properties of grain boundaries 
in G?-wave superconductors. We have demonstrated how 
they can support a superconducting state with broken 
time reversal symmetry. We established a connection be- 
tween such a state and the existence of vortices carrying 
a fractional flux. The key to this connection lies in the 
observation that a T-violating state at the grain bound- 
ary can lead to a non-trivial phase shift in the Josephson 
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current-phase relation. Vortices occur at locations where 
this phase shift changes. Corners of grain boundaries are 
important places for such vortices, because they separate 
segments with different properties. 

Our results compare qualitatively well with the exper- 
imental observation by Kirtley and coworkers p6|. At 
present it is, however, unclear whether we should con- 
sider these experimental results as an evidence for frac- 
tional vortices and, consequently, for the presence of a T- 
violating superconducting phase. Obviously, if the length 
scale of the magnetic field along the grain boundary is 
comparable with the distance between the corners or the 
vortices, then it is impossible to associate a definite flux 
with each vortex separately. In order to draw a firm con- 
clusion, this point has to be clarified experimentally. 

The conditions for the observation of this effect are 
best if the crystal orientations of the grains are chosen 
so that the grain boundary faces a lobe of the d x i_ y i- 
wave pair wavefunction on one side and nearly a node 
on the other. As we pointed out in Sect. Ill the lat- 
ter boundary provides a particularly good situation for 
pair breaking of the d x 2_ y 2-wave component, which is 
an important condition for our scenario. Similar conclu- 
sions where found with alternative mechanisms p8| , p9| . 
Unfortunately, grain boundaries of this type are intrinsi- 
cally difficult to produce as homogeneous interfaces and 
are often hampered by irregularities. 

In conclusion we would like to emphasize that the ef- 
fect discussed here is a consequence of the exotic nature 
of the superconducting order parameter. No analoguous 
effect is possible in the case of a conventional s-wave su- 
perconductor, because grain boundaries have little effect 
on this pairing type. 
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financial support from the NSF-MRSEC grant DMR-94- 
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